source('dados.R')

hist(Danos, breaks=seq(min(Danos), max(Danos), length=6),
	right=FALSE, freq=FALSE,
	xlab='Danos', ylab='Frequência', border=rgb(0,0,1,.65), lwd=2, col=rgb(0,0,1,.05),
	main=NULL)

X <- seq(0, 180, by=10)
lines(X, dlnorm(X, mean(log(Danos)), sd(log(Danos))), col='green', lwd=2)
lines(X, dgamma(X, (mean(Danos)^2)/var(Danos), mean(Danos)/var(Danos)), col='purple', lwd=2)

legend('topright', c('Log normal', 'Gama'), lwd=rep(2, 2), col=c('green', 'purple'), cex=.8)

